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ABSTRACT 

We present a [60-100] vs. [100-170] /itm two colour diagram for a sample of 61 young 
stellar objects (YSOs) observed with the Long Wavelength Spectrometer (LWS) on- 
board the Infrared Space Observatory (ISO). The sample consists of 17 Class sources, 
15 Class I, 9 Bright Class I (L bo i > 10 4 L & ), 20 Class II (14 Herbig Ae/Be stars and 6 
T Tauri stars). We find that each class occupies a well defined region in our diagram 
with colour temperatures increasing from Class to Class II. Therefore the [60-100] vs. 
[100-170] two colours diagram is a powerful and simple tool to derive from future (e.g. 
with the Herschel Space Observatory) photometric surveys the evolutionary status 
of YSOs. The advantage over other tools already developed is that photometry at 
other wavelengths is not required: three flux measurements are enough to derive the 
evolutionary status of a source. As an example we use the colours of the YSO IRAS 
18148—0440 to classify it as Class I. The main limitation of this work is the low 
spatial resolution of the LWS which, for some objects, causes a high uncertainty in 
the measured fluxes due to the background emission or to the source confusion inside 
the LWS beam. 

Key words: stars: circumstellar matter - stars: formation of - stars: fundamental 
parameters - stars: pre-main sequence 



1 INTRODUCTION 

It is generally accepted that the evolutionary status of young 
stellar objects (YSOs) is characterised by the shape of their 
spectral energy distribution (SED). During the early phases 
of star formation the emitted SED has the shape of a cold, 



few ten 3 of K, blackbody modified by the emissivity of the 
dust. The central object, even when already tormed, is not 
visible at any wavelength since the surrounding dust en- 
velope completely obscures the inner region. As the evolu- 
tion goes on, the optical depth of the dust decreases making 
the central object optically visible until its spectrum looks 
like that of a normal star with an infrared excess due to 
the residual circumstellar material. It is then clear that the 



* ISO is an ESA project with instruments funded by ESA Mem- 
ber States (especially the PI countries: France, Germany, the 
Netherlands and the United Kingdom) and with the participa- 
tion of ISAS and NASA. 



characterisation of the evolutionary status of YSOs is better 
done by observing the dust emission rather than the stellar 
radiation which is only visible during the late phases of star 
formation. 



ISO, the Infrared Space Observatory (Kessler et al 



1996 ), has opened a new observing window covering the near 
to far infrared wavelengths up to A ~ 200/im, where the dust 
emission peaks. It is then natural to use ISO data to derive 
the evolutionary status of YSOs. To this aim we introduce 
a two colour diagram defined between 60 and 170 /im which 
exploits almost all the spec tral coverage of t he Long Wave- 
length Spectrometer, LWS ( |Clegg et al. 1996] ) . Already with 
IRAS data two colour diagrams h ave been used to infer the 
evolutionary phase of YSOs, (e.g. Beichman et al. 1986, or 



Berrilli et al. 1989), however, the LWS extends the FIR (Far 
InfraRed) spectral coverage of IRAS by a factor ~ 2 in the 
region where the peaks of emission of blackbodies with tem- 
peratures in the range 15 — 30 K fall, typical of YSOs in the 
early phases of their evolution. 
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Tommasi & Pezzuto (1997) and Pezzuto et al. (199£ ) 



have indeed already shown that such a diagram allows to 
identify the evolutionary status of YSOs. While in the pre- 
vious work a small sample of objects was used, we have now 
completed the analysis of the whole sample of YSOs ob- 
served by the LWS as part of the core program dedicated to 
the study of the star formation. 

With respect to previ ous evolutionary diagnostic tools, 
like the spectral index a (Adams, Lada & Shu 1987), the 
bolometric temperature ( MyeTs^TT^adcT 1993|, Ghcn et al 
1995), ffld the Lboi vs. J\.3 mm diagram ( jSaraccno et af 



1996), our method can be applied to YSOs in all evolution- 
ary phases (Class 0, I and II objects), and does not require 
the knowledge of either Lboi, or the complete SED. 

In Section 2 we present the sample, and the data anal- 
ysis is briefly reported; in Section 3 the two colour diagram 
is shown and discussed; the conclusions are summarised in 
Section 4. 



2 OBSERVATIONS AND DATA REDUCTION 

Our sample consists of 61 YSOs for which an LWS full range 
grating spectrum, 43 — 197/xm with a resolution A/AA ~ 
200, was obtained. It corresponds to almost all the observa- 
tions c arried out with the LWS as part of the ISO-LWS core 



prograr ime on star formation, excluding a few sources whose 
spectra are too noisy to derive reliable photometric data. 
We have decided to include for this paper only the sample 
of our guaranteed time, to have a homogeneous set of data 
taken with the same observational procedure and compara- 
ble integration times. Our original sample, to which the 61 
objects discussed here belong, was selected on the basis of 
the knowledge of the instrument (spatial and spectral reso- 
lution, sensitivity) and on the spectral characteristics of the 
known sources, at the time of the definition of the observing 
programmes (around 1995). 

The raw data have been processed with version 8 of 
the off-line pipeline^ and analysed with the ISAP package. 
After removing bad points, the scans were averaged. The 
strong source correction, which takes into account the non- 



linearity of the detectors at high fluxes ( Leeks et al. 1999 
Leeks ioOOj ), was applied to the data if it was necessary. 



The LWS spectra are affected by instrumental effects 
(fringes, memory effects, bad flat field and dark current re- 
moval, see Gry et al. 2001) which might potentially influence 



the shape of the spectrum. Such effects can be, at least par- 
tially, corrected by using specific routines (LIA, the LWS 
Inter active Analysis). However it has been verified ( Grillo 
2001) that for the majority of our objects, the fluxes are not 
considerably changed by applying these routines. In the end, 
we adopted the fluxes computed without any correction. 

For some objects, one or more off-source positions have 
been observed to estimate the bac kground emission. In our 
previous work (Pezzuto et al. 1999) we compared the colours 
before and after removing the background and we found 



1 After submitting the manuscript, the final version of this soft- 
ware has been released (version 10). The changes in the data pro- 
cessing with respect to version 8 are not relevant to the present 
work, as we verified for some faint or noisy sources (for instance 
SR. 9). 



that, again in a statistical sense, the subtraction does not 
alter our results. Moreover, in star forming regions the pres- 
ence of a diffuse, not homogeneous, cloud emission makes it 
difficult to define the background level. In addition, when 
the instrument does not have enough spatial resolution, the 
confusion due to other sources in the beam can make even 
the meaning of background dubious. For these reasons we 
did not subtract the background from the on-source spec- 
tra. The large spread of the colours in our diagram is, likely, 
the consequence of this choice, as we discuss in the next 
section. 

The colour [Ai — Aj] is defined as: 
[A 4 -A j ]=log J 

Fij being the flux densities in Jy, averaged over a certain 
band, at the wavelengths Xij. In the case of photometric 
data A is the effective wavelength of the adopted band. Deal- 
ing with spectra the choice of the wavelength is somewhat 
arbitrary. A natural choice is to measure the flux at 60 /im 
and 100 /im, which allows a comparison with IRAS data. As 
a third wavelength we adopted 170 /im: this choice was first 
motivated by t he instrumental characteris tics of the LWS 
(as discussed in Tommasi & Pezzuto 1997). It then turned 



out that these three wavelengths are close to the effective 
wavelengths of the broadband filters of PACS (Poglitsch & 



Waelkens 2001), a far infrared camera for imaging and spec- 



troscopy which is part of the payload of the future ESA cor- 
nerstone mission the Herschel Space Observatory. In such a 
way it will be possible to use our two colour diagram with 
PACS data. 

For the bandwidths there is no obvious choice. At first 
sight a good solution could be to simulate the spectral re- 
sponse of the IRAS bands, as the ISAP package allows. 
There are, however, at least two reasons for not following 
this approach. The first one is that this choice is not applica- 
ble to the third band since IRAS did not observe at 170 /im. 
Secondly, future space missions will enlarge the number of 
known YSOs, so it is important that our tool does not de- 
pend on the photometric bands adopted in previous observa- 
tions. So, we simply averaged the spectra over 1 /im around 
the three wavelengths, and verified that our conclusion are 
not altered when data derived with a completely different 
instrument, in particular with IRAS, are used. This topic 
will be addressed in the next section. 

The uncertainty of the LWS fluxes can be as low as 10% 
for on-axis, not too bright nor too faint, point sources. These 
conditions are seldom met for our sources, so we adopted a 
more conservative error of 30%, which turns into an error 
on each colour of 0.26. 



3 RESULTS AND DISCUSSION 

In Table ji] we report, for each object: its name; the identifier 
in the IRAS Point Source Catalogue, if present; the LWS 
and, if available, the IRAS colour [60-100]; the LWS [100- 
170] colour. The membership of each source to a class has 
been assessed by looking at the literature. 

The [60-100] vs. [100-170] diagram is shown in Figure 
For comparison we reported the colours, computed following 
the same procedure outlined before, of a blackbody at vari- 
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Table 1: The sample of YSOs observed. Sources have been grouped according to their known evolutionary status. Inside each group 
sources have been sorted in increasing right ascension. An asterisk flags objects classified as FU Orionis. The [60-100] IRAS colour has 
been computed by using the fluxes reported in the SIMBAD database or the fluxes given in the quoted reference. The IRAS colour of 
L 1448 N has been obtained by adding the fluxes of L 1448 N and L 1448 NW, which fall inside the LWS beam. 
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ISO-LWS colour-colour diagram 



1.0 




1 








V 










A 


Class I 


o 


- 


0.8 












- >K 


Bright Class I 




- 




• 


HAEBE stars 




- 


0.6 


— 




„ <P o 






O 


T Tauri stars 


yC^O 9 Q ^ o 

/ o 


o- 








0.4 


— 






— 


0.2 






x*. 4 *a a 4 
















3X 






0.0 




V o 2 


* • 
















>< * 






0.2 




, 







-0.4 -0.2 0.0 0.2 0.4 

[100-170] 

Figure 1. Two colour diagram for the sample of YSOs. Bright Class I sources are those with bolometric luminosity Lbol > 10 4 Lq. 
Boxed symbols identify FU Orionis objects. The error, not displayed, on each point is 0.26 in both colours, assuming an uncertainty on 
the fluxes of 30%. Solid line: blackbody colours at temperatures (marked by crosses): 25 K (upper right corner), 30, 40, 50, 60, 70, 80 
and 90 K (lower left corner). The numbers identify sources which are individually discussed in the text: 1) MWC 297; 2) FIR 5 and 3) 
FIR 3 (both in NGC 2024); 4) L 1206; 5) IRAS 2 (in NGC 1333); 6) NGC 281 W; 7) IRAS 18148-0440; 8) IRAS 13547-3944; 9) SR 9; 
10) T CrA; 11) IRS 44 and 12) IRS 43 (both in p Oph); 13) HH 26 IRS. 



can look deeper into the inner, and warmer, regions of the 
envelo pe and, consequently, T c increases. 



[100 — 170] colour and falls closer to Class I sources. It is 
report ed as a Class in the co mpilation by Andre et al 



Ths remaining of this section is organised as follows: 
first we discuss separately each class and argue that the 
presence of outlyers, evident in the diagram, can be ten- 
tatively explained, in many cases, as due to background 
contamination. For such reason we cannot average all the 
colours for each class because this procedure assumes that 
the dispersion of the points is caused only by random fluctu- 
actions. Instead, we preferred to compute the median which 
is less sensitive than the mean to the outlyers. Then, in Sec- 
tion 3.4, we test our result with a different data set, using 
IRAS fluxes, to verify that we did not introduce any bias in 
the data analysis, at least for the [60 — 100] colour. Finally, 
the medians are reported in Figure ^| along with their root 
mean square (RMS) errors, and their position with respect 
to the blackbody colours is discussed in Section 3.5. 



3.1 Class sources 

Class sources appear well clustered in the upper right 
region of the diagram with [60 — 100] > 0.45 and [100 — 
170] > 0.08, except four objects discussed below. The me- 
dian "colours are: [60 - 100] = 0.53 ± 0.23 and [100 - 170] = 
0.19 ±0.18. 

IRAS 18148-0440 (number 7 in Figure @) has a smaller 



(2000| ) but [Tafalla et al. (200C| ) argue that this source is 



already in a transition state between Class and I. From 
our diagram we suggest that IRAS 18148 would be better 
classified as Class I. 

The two sources FIR 3 and FIR 5 in NGC 2024 (num- 
bers 2 and 3 in the figure) have blue colours, similar to those 
of Bright Class I (Lbol > 10 4 £ a). From the analysis of the 
same spectra used in our work, |Giannini et al. (2000[ ) have 
shown that in NGC 2024 the emission lines are mainly pro- 
duced by the extended cloud in which the sources are embed- 
ded. Since the Hn regio n associated with N GC 2024 has a 
luminosity of 2.8- 10 4 L ( [Mezger et al. 1992| ), while for FIR 5 
Wiesemeyer et al. (1997) estimate Lboi ~ 1 — 41/0, from the 
position of the sources in our diagram we suggest that the 
observed continua are dominated by the diffuse PDR more 
than by the objects themselves. 

The source IRAS 2 in NGC 1333 (labelled 5 in the plot) 
falls close to Class I I objects. The IRAS colour found by 
Jennings et al. (1987) and reported in Table ^ is redder than 
ours and puts the source closer to other Class 0. The fluxes 
were measured with the Chopped Photometric Channel of 
IRAS which provided a higher angular resolution than that 
specific to the IRAS survey, but it is worth noting that the 
colour directly measured by IRAS is 0.35, not very different 
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from our value of 0.30. It is then likely that the colour we 
found is strongly affected by the low spatial resolution of 
the LWS rather than depending on intrinsic characteristics 
of the source. 



3.2 Class I and Bright Class I sources 

Without considering the two FU Ori sources, we find for 
Class I objects the following median colours: [60 — 100] = 
0.280±0.075 and [100-170] = 0.09±0.10. The large scatter 
in the [100 — 170] colour is due to the position of L 1206 
(source number 4) close to the Class II objects, and to a little 
"cluster" of three sources with [100 - 170] > 0.2: for IRS 43 
and 44 (numbers 12 and 1 1 in the figure) the large [100 — 1 

k 



Table 2. Comparison between median [60 — 100] colours derived 
with LWS and IRAS data. 



ZQ]_ 



colour ]s due to the diffuse emission in the p Oph complex 



in which they are embedded. In fact, the subtraction of the 
background moves their colours into the Class I region. 

Objects having a total luminosity L > 10 4 Lq have been 
considered as a distinct group (Bright Class I), since such a 
high luminosity is usually associated with very massive hot 
young stars whose FIR colours are expected to be different 
from those of lower luminosity Class I sources. In fact, with 
the exception of NGC 281 W (number 6), they have [60 — 
170] < 0.2 and [100 - 170] < -0.1, and fall close to the 
blackbody line with colour temperatures T > 50 K. The 
median colours are [60-100] = 0.02 ±0.16 and [100-170] = 
-0.130 ±0.083. 

The col ours of NGC 281W, whose total luminosity is 
2.4 ■ 10 4 L @ ( |Carpenter, Snell fc Schloerb 1990[ ), are the red- 



dest of this group and are in strong disagreement with those 
of the other objects. We do not have any explanation for 
such a position. 



3.3 Class II sources 

This group is divided into low mass stars, M < 2 Mq , called 
T Tauri after their prototype, and intermediate mass stars, 
2 M @ < M < 8 Mq, called Herbig Ae/Be stars (or HAEBE). 
This division reflects physical and observational differences 
between them, but in our diagram they occupy the same 
region so, to increase the statistical significance of our work, 
we will discuss them together. All the Class II objects, with 
the exception of MWC 297, IRAS 13547, T CrA and SR 
9, cluster in a well defined region of the plot. The median 
colours of all sources, except the FU Ori, are [60 — 100] = 
0.11 ± 0.14 and [100 - 170] = -0.03 ± 0.12, with a colour 
temperature T c ~ 50 K. 

For SR 9 (number 10 in the figure) the fluxes in the on- 
source position are equal to those in two observed off-source 
positions so that we conclude that the continuum emission 
of this T Tauri star is completely dominated by the intense 
diffuse background of the p Oph region, in which SR 9 is 
embedded, with colours in agreement with those of Class 
sources (see Table [l]). Its [100 — 170] colour is similar to that 
of IRS 43 and 44 which also lie in the same region (see the 
previous discussion). 

The anomalous position of the HAEBE star MWC 297, 
in the bottom left corn er and labelled w ith 1, is compati- 



ble with the finding of Drew et al. (1997) that it could be 
in a more evolved stage and already on the main sequence. 
In the opposite corner, IRAS 13547 (number 8) and T CrA 





LWS 


IRAS 


Class 


0.53 ±0.23 


0.53 ±0.16 


Class I 


0.280±0.075 


0.385±0.078 


Bright Class I 


0.02 ±0.16 


0.13 ±0.16 


Class II 


0.11 ±0.14 


0.12 ±0.20 



(10) have colours compatible with Class sources: even if 
their membership to HAEBE is doubtful, the former be- 
ing defined as "potential candidate" and the latter reported 
among "other early type emission line stars with IR excess" 



in the compilation of HAEBE by The, de Winter fc Perez 



(1994), these sources are optically visible and they cannot 
belong to the group of Class 0. The only conclusion that we 
can draw from our diagram is that the FIR emission is either 
produce d by a diffuse backg round, for T CrA see the FIR 
maps of Wilking et al. (1985), or by a close cold companion. 



3.4 Comparison with IRAS data 

As we already pointed out, it is important that the result of 
our work, i.e. that the [60—100] vs. [100—170] diagram is able 
to separate YSOs belonging to different evolutionary stages, 
is independent of the instrument used and on the adopted 
photometric bands. Further, since we neglected instrumental 
effects, did not take into account the background contamina- 
tion, defined two colours chosing three arbitrary wavelengths 
and averaging the fluxes over three arbitrary bandwidths, 
the reader could wonder about the reproducibility of our 
finding with different instrument /samp le/data analysis. 

To address some of these points we compared the me- 
dian of the [60 — 100] colours we found, with the medians 
derived using the IRAS fluxes, when available in literature. 
These two sets of data have in common only the effective 
wavelengths: bandwidths and data acquisition/reduction are 
completely different. As in our original sample, we did not 
consider FU Orionis objects. The two sets of medians are 
reported in Table |^. 

For Bright Class I the number of sources with IRAS 
data is too small to draw any conclusion. For Class and 
II the agreement between the two sets of data is excellent. 
Only for Class I is the difference larger than 1 RMS. From 
our point of view, however, the difference between 0.280 and 
0.385 does not make a big difference in the conclusion that 
Class I sources occupy a specific region of the diagram. What 
this difference means is just that we cannot specify with 
great accuracy what the median [60 — 100] colour is. 



3.5 The median colours: dust emissivity vs. 
source confusion 

The median colours are reported in Figure |2| In this fig- 
ure we also plotted, as in Figure the colours of a single 
blackbody at different temperatures (solid line). In order to 
clarify why the YSOs have colours under this curve we also 
plotted the colours of a blackbody modified by an emissivity 
law oc A - * 3 which better approximates the emission proper- 
ties of the dust, the main emitter in the FIR. The dashed 
line corresponds to the case [3 = 2, typical of interstellar 
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Figure 2. The median colours for each class (filled symbols) with their root mean square errors (dashed lines). From top to bottom: 
Class 0, Class I, Class II and Bright Class I. Solid line: colours of a single blackbody as in Figure 1; dotted line: colours of a blackbody 
modified with an emissivity law oc A -1 at T: 25, 30, 40, 50 K; dashed line: as before but with the emissivity oc A -2 and T: 25, 30 K. 
Dotted-dashed line: colours of two blackbodies at 25 and 90 K and different radii; the cross marks the colours corresponding to the case 
when the radius of the colder source is 10 times larger than the radius of the warmer source. 



dust, while the dotted line f3 = 1, is more representative of 
the dust in the star forming regions. It is clear from the fig- 
ure, that no positive j3 can account for the positions of the 
median values. 

Another possible explanation is related to the confu- 
sion of the sources inside the LWS beam. It is not possible 
to model this effect since we should know, object by object, 
the physical properties and the spatial distribution of the 
sources around the observed position. However, we can at 
least check if this hypothesis moves the colours in the region 
of the plot where their medians lie, by simply summing the 
fluxes of two blackbodies. In this case there are three free pa- 
rameters: the two temperatures and the ratio of their radii. 
As an example, we have plotted in Figure |^, the dotted- 
dashed line, the colours obtained for the temperatures 25 
and 90 K, with a ratio varying from to oo (the cross sym- 
bol marks the case corresponding to the radius of the colder 
body 10 times larger than the radius of the warmer one). 
As can be seen the sum moves the colours in the wanted 
direction. 

The effects of a finite optical depth of the circumstel- 
lar envelope, while they cannot be accounted for by a sum 
of blackbodies, should however be in qualitative agreement 
with the previous model, so that we can expect also for iso- 
lated objects a displacement from the blackbody colours as 
observed in our plot. 



4 CONCLUSIONS 

The [60 — 100] vs. [100 — 170] two colour diagram has been 
presented for a sample of YSOs observed with the LWS in- 
strument on-board ISO. Each class of YSOs occupies differ- 
ent regions of the diagram with bluer colours corresponding 
to older classes. The median values of the colours are: 



Class 0: 



Class I: 



[60 - 100] = 0.53 ± 0.23 
[100-170] =0.19 ±0.18 

[60 - 100] = 0.280 ± 0.075 
[100 - 170] =0.09 ±0.10 



Bright Class I: 



[60 - 100] = 
[100 - 170] 



0.02 ±0.16 
= -0.130 ±0.083 



Class II: [60 - 
[100 



100] =0.10 ±0.16 
-170] = -0.04 ±0.12 



The relative position of the different classes is related to the 
physical properties of the circumstellar matter and is not 
an artefact of our data analysis. This has been tested, for 
the [60-100] colour, through a comparison with IRAS data 
which also show this segregation effect. It is not possible 
to make such test for the other colour. However it is very 
unlikely that the class separation is due to a bias either in 
the sample selection or in the data analysis. The change of 
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the colours is a genuine evolutionary effect: in the younger 
objects we can only observe the outer, and colder, regions of 
the dusty envelope. The median colours are systematically 
under the curve of the blackbody colours, and the values 
we found are probably biased by the combination of both 
the background contamination and the confusion of the ob- 
served sources in the LWS beam. The low spatial resolution 
of the instrument is indeed the main limitation of our result. 

In this respect this work will be greatly improved by fu- 
ture space missions which will observe with a higher spatial 
resolution. For instance, the multiband photometric surveys 
of the FIR camera PACS at ~60, 100 and 170 /im, will be 
able to make the census of each different evolutionary class 
of protostellar and pre-main sequence objects in the most 
important star forming molecular clouds. 

In spite of the present limitation, we conclude that this 
two colour diagram offers a powerful and simple tool for 
the evolutionary classification of YSOs. The advantage over 
other methods is that follow up at other wavelengths to iden- 
tify the evolutionary status will not be needed. As an exam- 
ple, even if we stress that with the LWS data our diagram 
works mainly as a statistical tool, with only three fluxes, 
and no other photometric measurement, it is possible to as- 
sign the evolutionary class to the source IRAS 18148—0440: 
Class I, in agreement with Tafalla et al. (2000), rather than 



Class source ( Andre et al. 200C ) 
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